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RESEARCH SUMMARY 


This paper presents formulae tables, and figures that canbe used to estimate 
the ratio of mean windspeed acting on the flame from a spreading wildland fire 
to the mean windspeed 20 ft (6 m) above the vegetation cover. The formulae are 
based upon the logarithmic windspeed variation law characteristic of constant 
shear turbulent flow, and are restricted to flat terrain with uniform, continuous 
vegetation cover. 

By equating integrated bulk drag due to crown foliage to the shear stress at 
the top of the crown layer, a model for relating windspeed within and below a 
uniform forest canopy to windspeed 20 ft (6 m) over the canopy is developed. 
Important variables in this model include stand height, crown closure, foliar 
surface/volume ratio, and crown bulk density. Comparison of model predictions 
with reported experimental measurements shows good agreement. 
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INTRODUCTION 


The behavior of fire in wildland fuels is influenced strongly by the wind acting 
on the fire. Predicting fire behavior for fire control planning or preparing burning 
prescriptions must include this influence. Rothermel (1972) gives a mathematical model 
for predicting the rate of spread of a surface fire in wildland fuels. This model uses 
an average windspeed "at midflame height'' to account for the influence of wind on the 
rate of spread. But windspeed is usually measured or forecast at a standard height 20 
ft (6 m) above the vegetation (Fischer and Hardy 1972), and varies with height and veg- 
etation cover. The poorly defined "midflame'' windspeed can be approximated by using a 
Spatially averaged value of the windspeed over an appropriate height range. 


Presented here is a means of estimating the windspeed over various wildland fuel 
types from the ''standard" windspeed at 20 ft above the vegetative cover. The report 
is divided into two major sections. The first deals with the wind above a vegetative 
cover that is a single-stratum fuel (grass, brush, etc.). The second part deals with 
wind under a forest canopy. The windspeed at the top of the forest canopy is found by 
use of a logarithmic wind profile. Windspeed beneath the canopy layer is estimated from 
the canopy-top value by use of a model based on mechanical force balance. 


This research report deals only with the steady, undisturbed windfield and its in- 
fluence on fire in surface fuels. No account is taken of the influence of the fire on 
the windspeed, wind direction, or the profile of windspeed with height. Furthermore, 
flat terrain and uniform continuous vegetation cover are assumed. 


iheScymestrictions clearly lami the applicability of the results, but even the 
extent and severity of the restrictions is uncertain. Rothermel's model implicitly 
accounts for local (near-flame) influence of the fire on the windfield. At what inten- 
Sity or size the fire's influence on the prevailing windfield becomes significant can 
only be estimated by order-of-magnitude arguments at present. Future research should 
help to clarify this situation. 


iheasrestriction to flat terrain can probably be relaxed to ''smooth"” or "slightly 
uneven" terrain, but numerical descriptions of adequate ''flatness'' are not presently 
available. Similar considerations apply to the uniformity and continuity of vegetation 
cover required, and these matters are under study in the field of applied meteorology 
(Bergen 1976). 


With these restrictions in mind, this effort is seen as a small first step in the 


direction of a more complete description of the complex phenomena with which we are 
dealing. 


LOGARITHMIC WIND PROFILE 


The windspeed above a vegetative cover was determined by using the logarithmic wind 
profile in the following form (Sutton 1953, p. 239): 


where 
U_ is the average windspeed at height z 


is, the friction velocity, (U, = ¥t/p, tT as the horizontal 
shear stress and op is air density) 


K = 0.4 (the von Karman constant) 
Zz is height above ground 
D_ is the zero-plane displacement 


Z5 is the roughness length. 


The values for Dy and Zz, were taken from the works of Cowan (1968) and Stanhill 
(1969). Figure 1, taken from Stanhill, shows a plot of Dy versus vegetation height. 
The plot agrees with the relationship Dy = 0.64H derived by Cowan, where H is the vege- 


tation height. The derived relationship was based on 2 0.13H as found by Tanner 
and Pelton (1960). 
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Fig. 1.--Relattionshtp between zero-plane displacement Dg and vegeta- 
tion height H for different vegetation types (after Stanhtll 1969). 


WINDSPEED AT THE TOP OF THE VEGETATION COVER 


To predict wildland fire behavior, the "standard" windspeed at 20 ft must be 
converted to a value representing the actual motion of the air in the vicinity of the 
fuel and flame structure. As a first step, we related the 20-ft windspeed (U5o44) 
to the windspeed at the top of the vegetation (Uy) - ri 


The relationship between U., and U5 94H 1S) chen: 
= P 20+H - 0.64H f ; 3 

Usoy =(U,/K) 2n =o oRET eas for H measured in feet. (2) 

Since 
= Q H — ee 0.644 = z = 

U, = (U,/K) an 0138 1.02 U.//K U,/K (3) 
we can write equation (2) as follows: 

meee cle (4) 

Wane ,,(20+0. 36H 


WINDSPEED FOR PREDICTING FIRE SPREAD 
IN SURFACE FUEL (NO CANOPY CASE) 


With the profile of windspeed variation defined, we need only choose an appropriate 
height range over which to average to obtain a relationship between 'midflame!'' windspeed 
and windspeed at 20 ft above the fuel surface. For this purpose here we assume that the 
flame extends to some height above the surface vegetation and average the windspeed be- 
tween the top of the vegetation and this height. 


On 


Consider an average windspeed over the height range H to H + He, where Hy. is the 
extension of the flame above the fuel bed. 


= er atl in ih as 2-D, 
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Note that at the lower limit (H) the value is nearly zero since: 
H-0.64H \ + 
in (4558 oe 


Thus we find 


We ilar dike Sy Os tn H + He - 0.64H 
U = E gn F =| 
K He 0.13H 
U, r H-/H + 0.36 
= aa ( + 0.36 7 on ARR OCR -] (8) 


Values of U can be calculated for various windspeeds at (20+H) ft over single stratum 
fuels. Using equation (2), we find: 


TT 1 + 0.36H/H ie lal ae Oe Ko 
== = oe eeee F Qn all gal eee =| (9) 
Wa 2 peer ZOLeNONSEHN | 0.73 
ORISh 


The graph of figure 2 can be used to establish the ratio of the 'midflame'' windspeed to 
the windspeed 20 ft over the vegetation cover for various fuel heights H and flame 
extensions H,. 
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Fig. 2.--Average windspeed acting on a flame extending above a uniform surface fuelbed 
layer (vegetation cover), due to log windspeed variation. 


GENERAL USE OF THE LOGARITHMIC WIND PROFILE EQUATION 


According to several authors, the log-wind profile represents the actual wind pro- 
file near the surface of the earth only under neutral stability condition. However, 
various researchers have found that the log-wind profile represents the wind above a 
vegetative cover (very rough surface) over a wide range of atmospheric conditions. For 
example, van Hylckama (1970) reports that a log-wind profile represented the wind vari- 
ation above a salt cedar stand (Zamarix pentandra) for all his afternoon observations 
and for many night conditions. Of 744 data sets, only 84 were too irregular for log- 
wind representation. Oliver (1971) found the wind profile above a pine forest to be 
well described by the logarithmic form over a wide range of stability conditions. 


Apparently, large temperature gradients cannot be maintained above a forest canopy 
due to air movement through the canopy (Munn 1966, p. 158) and due to mixing by which 
the wind blowing over a very rough surface reduces air temperature variations above the 
canopy. Thus extreme lapse or inversion temperature conditions are unlikely to obtain 
over a vegetative cover, resulting in a broad applicability of the log-wind profile. 
However, we advise that this form not be used to estimate winds under conditions of 
strong surface air temperature inversions. In summer this happens most often at night 
under clear skies with very low surface wind speeds. So when wind represents an impor- 
tant influence on fire behavior, the log profile usually is applicable. 


WIND UNDER A CANOPY 


To model the variation of windspeed with height for air flow through and under a 
forest canopy, we make several simplifying assumptions: 


1. Below some height, near but below the top (height H) of the uniform 


z 

M’ 
forest canopy, the windspeed is approximately constant with height. 

2. The foliage within the height range from Zu down to the bottom of the 
live crowns at height a provides a bulk drag force that resists the 


airflow. 


3. Shear stress (equal to that in the constant stress layer above) on the 
surface z = z,, balances the integrated bulk drag force in the constant 
speed layer. 


The approximate windspeed profile is thus: 


z - 0.64H 
(U,/K) 2n (2735 ) ; Z>Zy 


(10) 


Ue i ea 


The assumption of a constant windspeed beneath the canopy seems quite reasonable accord- 
ing to various authors (Fons 1940, Geiger 1966, Shaw 1977). Now the resistance to 
airflow in the canopy can be expressed as a force/unit volume, R, due to drag on the 
canopy foliage (conifer needles in most applications). The value of R is given by: 


ei a 


where 
0 is the air density 
U_ windspeed in canopy 
C. the average drag coefficient (random foliage orientation) 
A the side view area of needle = 2D (length x diameter) 


n is the average number of needles per unit volume 


n ee = 8B, the average "packing ratio" (Rothermel 1972) (12) 
thus 
, = 
nA =nT r? : = = ee (13) 


where g is the surface-area-to-volume ratio of a single needle. This bulk resistive 
force, integrated over the height range Zz to 2, provides the shear stress requirement 
(1) for the upper surface: 


t=R “(Zy - Zam) (14) 
Since at Zu the shear stress is assumed to be that of the constant-shear layer, we have 
t= 9 US (15) 
and so 
(2 = za)? Uc Cr Bo/27 = p U2 (16) 


From this last expression and equation (3) we find 


(u,/u,) = 2nK2/ (, Bo (2, - 7) (17) 


For a uniform stand of nearly identical trees we can replace 2y 72 by H-(CR) where H 


is tree height and CR is the crown ratio. Likewise, the average packing ratio, 8, 
can be expressed as 


B= ke (18) 


where 8 is the packing ratio for a single tree crown and F is the fraction of the canopy 
layer filled with tree crowns. The fraction F is to be approximated by the product of 
crown closure and a fraction accounting for the tapering of crowns that results in addi- 
tional void volume higher in the canopy. Finally, the product of F and CR can be 
represented by f, a fraction that represents the portion of volume under the canopy 

top that is filled with tree crowns. 


Typical crown ratios! are displayed in table 1. The fraction of the canopy layer 


that is occupied by tree crowns (F above) we take to be 0.4 for dense stands and 

0.1 for open stands. Table 2 gives the volume-filling fractions (f above) resulting 
from the products of the crown ratios of table 1 and these closure factors. In terms 

of the volume fraction, f, equation (17) provides a relationship for the windspeed under 


the canopy, Ue» and the windspeed at the top of the canopy, Un: 


U/ Uy =< (2n/tsotn) (19) 


Table 1.--Typtcal crown ratios, percent 


Stand : Shade-tolerant trees : Shade-intolerant trees 
stocking : Young : Mature : Young : Mature 
Dense 80 60 40 20 
Open 90 70 70 50 


Table 2.--Volwne filling fractions (factor f), percent 


Stand : Tolerant : Intolerant 
stocking : Young : Mature : Young : Mature 
Dense 32 24 16 8 
Open 9 7 7 5 


This expression can be reduced to a simpler form by using nominal numerical values for 


the parameters that are essentially constant. The value of K can be taken as 0.4, Cy abs 


approximately 1.0 (Schlichting 1968), and the product go can be treated as a constant? 


equal to 3.26 Et. So for H measured in feet, 


oe ee eee. 0.556 on 
Uy Vo Bo fH (1) (3.26) fH JfH 


Values of U /Uy for three tree heights are given in table 3, based on the f values of 
table 2. 


1 The authors are indebted to Dr. James K. Brown of the Northern Forest Fire Labora- 
tory for providing these typical quantities. 

2 This is an average value for 10 Rocky Mountain conifers. Data were taken from study 
material by J. K. Brown (on file at NFFL, Missoula) used in establishing crown foliage 
weight relationships (see Brown 1976). The standard deviation of the product 80 for 
these species is 1.08 fen 


Table 3.--Windspeed ratto U Uy 


: Tolerant species : Intolerant species 
Stockungay: Young wy MatuxeyeyPacatich: Young : Mature : Pacific 
level : AO OOM: SOM ie: AQ) BLOOM See AUSKOY 
Dense 0.155 (Oe, UAL 0.084 OF2205) (O96 0.146 
Open 0.293 0.210 OUS6 ORSS2 0.248 0.185 


Three "typical" tree heights are used in this table to represent 
young stands (40 ft), mature stands (100 ft) and old-growth Pacific 
Coast stands (180 ft). For the latter, the crown ratios of "mature"! 
stands from table 2 were used. 


Since U. applies (almost) all the way to the ground, it is the 'midflame" windspeed. 


Hence the "midflame" windspeed in relation to the windspeed 20 ft over the canopy can 


be found by using equation (4) and choosing the ratio U SU, from table 3 or calculating 
it from equation (20). 


The equation for calculating UW. for stands with arbitrary values of f and H is: 


U/Uooiy = sf en (20 + 0.36H) 10.138) (21) 


The ratio UU 56414 is plotted in figure 3 for the typical and extreme values of f 
(Ealbilie=s2)) : 
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Fig. 3. Ratio of windspeed within (and below) forest canopy to windspeed 20 ft above 
canopy top. 
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COMPARISON WITH MEASUREMENTS 


Table 4 compares data reported by Shaw (1977) and Berglund and Barney (1977) with 
our calculated windspeed ratios U/U,, from equation (20) using the values of f from 


table 2. Table 5 compares calculated values of U/U from equation (21) and table 2 


20+H 
with data taken from Fons (1940), Raynor (1971), and Allen (1968). The agreement be- 
tween model results and measurements is surprisingly good in all cases. 


Table 4.--Windspeed ratto UL Uy 


WU 
: : : , os a from 
Species : Stand description : Data source : Calculated : published data 
Sitka spruce H = $455 £t, S.t., dense Shaw Gy77) 0.167 ORIG 
Scots pine nl = OSs) see5 Soils, Oeil Shaw (1977) 0.259 0.24 
Ponderosa pine H = 70 Ie Holla, Coe Shaw GS) ORZ97/ 0.34 
Japanese lanchweh—e SAMs te Se elemOp en Shaw (1977) 0.359 0.36 
Black spruce Il = Iss ste, “Sadlo,, OjoGM Berglund & 0.619 OFOS 5 


Barney (1977) 
*p. 10, Table 2, ratio of slopes of regression equations relating windspeeds within 
and above the stand. 


Table 5.--Windspeed ratto UJ Vso sH 


U_/U 


c. 20+H 
. . . . from 
Species : Stand description : Data source : Calculated : published 
i : : : data 
Ponderosa pine TO, its Selle 5 OEM Fons (1940) 0.185 ave. 0.182 
Red & white pine 34.5 ft, S.I., dense Raynor (1971) OF IBKS ave. 0.145 
Japanese larch Ns Ati, Solis, Gocom Alien (ISS) 0.180 ave. 0.147 
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APPLICABILITY OF RESULTS 


The numerical results derived in this report are subject to the restrictions listed 
below. The reader is cautioned to be certain that the specified conditions prevail 
when using these results. 


1. Flat terrain has been assumed throughout. If the terrain has substantial 
slope or roughness, the windfield will reflect this fact, and substantial deviations 
from these results may occur. 


2. Adequate fetch to establish a uniform friction layer has been assumed. The 
definition of ''adequate'' fetch is the subject of current research (Shaw 1977), so 
numerical limits cannot be stated at this time. But near forest edges, lakeshores, 
or transitions in surface vegetation cover, the results given here may not be accurate. 


3. A "well-behaved" windfield is modeled by the relationships herein. If the 
windfield is not steady, but fluctuates significantly in speed, direction, or both, 
the friction layer will be continually in a transient state as it responds to the 
forces at play. During the periods of such variability the results given here may not 
be applicable. 


4. Any interaction between the fire and the windfield that substantially influ- 
ences the speed or direction of the wind should invalidate these results. We have 
dealt here with a windfield whose structure is dominated by the influence of the vegeta- 
tion cover on the friction layer and any factor that disturbs this condition negates 
the validity of the results given here. 


Wal 
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